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We measure the branching fraction of B~ — > r~v-r using the full T(4S) data sample containing 
772 x 10 6 BB pairs collected with the Belle detector at the KEKB asymmetric-energy e + e collider. 
Events with BB pairs are tagged by reconstructing one of the B mesons decaying into hadronic final 
states, and B~ — > r~v T candidates are detected in the recoil. We find evidence of B~ — > t~v t with 
a branching fraction of B(B~ — > r~v T ) — [0.72lg '^(stat) ±0.11(syst)] x 10 -4 , and a significance of 
3.0 standard deviations including systematic uncertainties. 

PACS numbers: 13. 20. He, 14.40. Nd 



The purely leptonic decay B~ — > r~v r |l| is of par- 
ticular interest since it provides a unique opportunity to 
test the Standard Model (SM) and search for new physics 
beyond the SM. A recent SM estimation of the branching 
fraction based on a CKM fit is (0.733tg:£?3) x 10" 4 @. 
In the absence of new physics, a measurement of B~ — > 
t~v t can provide a direct experimental determination 
of the product of the B meson decay constant and the 
magnitude of the Cabibbo-Kobayashi-Maskawa matrix 
element /s|T4b| 01 ■ Physics beyond the SM, however, 
could significantly suppress or enhance B(B~ — > t~v t ) 
via exchange of a new charged particle such as a charged 
Higgs boson from supersymmetry or two-Higgs doublet 
models HQ. 

Experimentally, it is challenging to detect the B~ — > 
t~v t decay because it involves more than one neutrino 
in the final state, and therefore cannot be kinematically 
constrained. At e + e~ B factories, one can reconstruct 
one of the B mesons in the e + e~ T(4S) BB reac- 
tion, referred to hereafter as the tag side (Stag), either in 
hadronic decays or scmilcptonic decays. One then com- 
pares properties of the remaining particle(s), referred to 
as the signal side (B s i g ), to those expected for signal and 
background. The method allows us to suppress strongly 
the combinatorial background from both BB and contin- 
uum processes. 

The Belle collaboration reported the first evidence of 
B~ — > t~v t with a measured branching fraction of 
[1.79±g;^(stat)±g;f?(syst)] x 1CT 4 , and a significance of 
3.5 standard deviations including systematic uncertainty 
(a), using hadronic tags on a 449 x 10 6 BB data sam- 
ple @ . This was followed by measurements by Belle us- 
ing the semileptonic tagging method Q , and also by the 



BaBar collaboration using both hadronic [8j and semilep- 
tonic .3] tagging methods. The four results are consis- 
tent, and their avera ge b ranching fraction is found to be 
(1.67 ± 0.30) x 10" 4 [10j], which is higher than the above 
SM prediction based on the CKM fit, although the sig- 
nificance of the deviation is less than 3cr. Therefore, it is 
of great importance to improve the accuracy of the mea- 
surement, by increasing the luminosity and the detection 
efficiency. 

In this paper, we present a new measurement of B~ — > 
t~v t using a hadronic tagging method and the full data 
sample of the Belle experiment. The analysis described 
here has a number of significant improvements. The ma- 
jor improvements are from an increased data sample (a 
factor of 1-7), significantly improved hadronic tagging ef- 
ficiency (a factor of 2.2), and improved signal efficiency 
due to less restrictive selection requirements (a factor of 
1.8). The combined effect of these improvements and the 
accompanying change in the signal to background ratio 
due to the looser selection criteria results in a reduction 
of the expected error by a factor of two. The new analysis 
has also improved systematic uncertainties. 

We use a 711 fb _1 data sample containing 772 x 10 6 BB 
pairs collected with the Belle detector [ljj at the KEKB 
e + e~ collider operating at the T(4S') resonance [LH ]. 
Most of the data sample has been reprocessed using 
improved track finding and photon reconstruction. We 
use a dedicated Monte Carlo (MC) simulation based on 
GEANT [l3| to determine the signal selection efficiency 
and study the background. In order to reproduce the 
effect of beam background, data taken with random trig- 
gers for each run period are overlaid on simulated events. 
The B~ — > t~v t signal MC events are generated by the 
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EVTGEN package 14| with the radiative effects based 
on the PHOTOS code jl5|. To model the backgrounds 
from e + e~ — > BB and e + e~ — > qq (q = u, d, s, c) pro- 
cesses, we use large MC samples corresponding to about 
10 times and 6 times the luminosity of the data sample, 
respectively. We also use MC samples of rare B decay 
processes containing charmless hadronic and leptonic de- 
cays and b — s> u semileptonic decays, where the sizes are 
expected to be about 50 times and 20 times the data 
sample, respectively. 

The Bt a g candidates are reconstructed in 615 exclu- 
sive charged B meson decay channels using an improved 
full-reconstruction algorithm [lj|. The output full re- 
construction quality variable A/"t a g ranges from zero for 
combinatorial background and continuum events to unity 
if an unambiguous Stag is obtained from the hierarchi- 
cal neural network. We also use the energy difference 
AE = EB tas — Ecm/% and the beam-energy-constrained 
mass Mbc 



x10 J 



\PB t: 



where Eqm is the 



e + e~ center-of-mass (CM) energy, i?B tag an d PB tag are 
the energy and the momentum, respectively, of the -B tag 
candidate defined in the CM frame. Charged B tag candi- 
dates with Mag > 0.03, -0.08 GeV < AE < 0.06 GeV, 
and 5.27 GeV/c 2 < M hc < 5.29 GeV/c 2 are selected. 
Figure [T] shows a comparison of the M^ c distributions 
for different data sets and tagging methods. The tag ef- 
ficiency (0.24%) and the purity (65%) are improved by 
factors of 1.7 and 1.2, respectively, compared to Ref. [6J. 
The number of Stag's obtained for the full data set is 
1.8 x 10 6 . Note that in case of the B~ — > t~D t signal, 
in which the BB event has lower than average particle 
multiplicity, the tag efficiency is 2.2 times higher than in 
Ref. and is 0.31%. 

In events where -B tag candidates are reconstructed, we 
search for B~ — > t~v t decays. The r~ lcpton is iden- 



tified in the 



and 



de- 



e V e V T , (J, V^V T , 7T 

cay channels. Candidate events arc required to have 
one track with charge opposite to the -Btag candidate. 
The charged tracks are required to satisfy dz < 3 cm 
and dr < 0.5 cm, where dz and dr are unsigned im- 
pact parameters relative to the interaction point along 
the beam axis and the perpendicular plane, respectively. 
Charged tracks are classified as electron, muon, and 
pion candidates after rejecting kaon and proton candi- 
dates [ll|. Candidate r~ — > it'i^v-,- events are required 
to have one 7r° candidate reconstructed from 7r° — > 77 
for which no daughter photons are used in the -B tag 
reconstruction. The invariant mass of the ir~ir° state 
is required to be within 0.15 GeV of the nominal p~ 



mass [17| . Multiple neutrinos in the final states are 
distinguished using the missing mass squared variable 
M miss = ( e cm - E Btas - E Bsi J 2 - |p BtaB +pb 3 J 2 , where 
Eb s1s and PB sig are the energy and the momentum, re- 
spectively, of the B s [ g candidate in the CM frame. To 



avoid potential backgrounds from e 
and p~K L , we require M^ ias > 0.7 (GeV/c 
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FIG. 1: Mbc distributions for the -Stag candidate events. The 
triangles, open circles, and solid circles represent the distri- 
butions obtained by applying the original tagging algorithm 
in Ref. [f| to the previous data set, an improved tag to the 
previous data set, and an improved tag to the full reprocessed 
data set, respectively. The solid and dotted curves show the 
sum and the background component, respectively, in a fit to 
the full data sample. 



After removing the particles from the B tag candidate 
and the charged tracks and 7r°'s from the B s i g candi- 
date, there should be no other detected particles. We 
require that there are neither extra charged tracks with 
dz < 75 cm and dr < 15 cm nor extra ir° candidates 
("tt° veto") nor K\ candidates veto") 



llj. We 



define the extra energy -Eecl, which is the sum of the 
energies of neutral clusters detected in the electromag- 
netic calorimeter that are not associated with either the 
Stag or the n° candidate from the t~ — > it~it v t decay. 
The signal should have either zero or a small value of 
Eecl while background events have larger values due to 
the contributions from additional neutral clusters. For 
neutral clusters contributing to -Becl, we require a min- 
imum energy threshold of 50 MeV for the barrel and 100 
(150) MeV for the forward (backward) end-cap calorime- 
ter. A higher threshold is used for the end-cap calorime- 
ter because the effect of beam background is more severe. 
We require £"ecl < 1-2 GeV as the correlation between 
-Becl and M^ iss is small for each background compo- 
nent in this region. The B s i g reconstruction efficiency is 
improved by a factor of 1.8, compared to Ref. Q. This 
improvement mainly comes from looser requirements on 
the track momentum, and from an improved ir° veto, in 
which photons used in the -B tag reconstruction are not 
used in the veto. 

The signal detection efficiency is estimated based on 
MC samples after applying a correction for the -B tag re- 
construction efficiency. The correction factor is obtained 
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by fitting the Mb c distribution for an Eecl sideband 
sample denned by 0.4 GcV < Eecl < 1-2 GeV, for which 
the kinematics is expected to be similar to the signal. 
The resulting efficiencies are summarized in Table HI The 
validity of the efficiency estimation is checked by using a 
semilcptonic decay sample in which B s i g is reconstructed 
in the decay chain B~ — > D*°£~Dg (£ = e or /i) followed 
by £>*° -> D°tt° and D° R-7T+. 

The signal yield is extracted from a two-dimensional 
extended maximum likelihood fit to Eecl and M^^. 
The likelihood is 

t=i j 

where j is an index for the signal and background con- 
tributions, rij and fj are the yield and the probability 
density function (PDF), respectively, of the jth com- 
ponent, Ei and Mf are the Eecl and M^ iiss values in 
the ith event, respectively, and N is the total number of 
events in the data. We construct PDFs by taking prod- 
ucts of one-dimensional histograms in Eecl and M^ iss 
obtained from MC, except for the contribution of cross- 
feed from other B~ — > t~v t decays in B~ — > t~v t , 
t~ — > tt~v t . The signal component in t~ — > ■k~v t 
candidate events includes large cross-feed contributions 
from — > l~viv T and r~ — > ir~ir°v T decays. A two- 
dimensional histogram PDF is used for the latter con- 
tribution to take into account the correlation originating 
from the misreconstructed ir°. The dominant background 
components are from BB decays to a final state with 
charm and continuum processes. The small background 
from charmless B decays is also included in the fit. In 
the final sample with Eecl < 1.2 GcV, the fractions of 
the background from B to charm decays, continuum, and 
charmless B decays are estimated from MC to be 89.8%, 
0.5% and 9.7% for leptonic r decays and 75.1%, 18.4% 
and 6.5% for hadronic r decays. 

Other B decays from which only one charged parti- 
cle is detected can make a peak near zero Eecl and 
mimic the signal. The dominant peaking decay modes 
arc B~ -)• D^ £-D t and B° -)• D^+i've, decays 
where the D decays scmilcptonically or to a final state 
with one or more K^'s. Charmless B decays such as 
B~ — > 7i £~i>i, K~vD, KjjK~ , -ftT*~7, and ^"^7 also 
contribute. The total fraction of these B decay modes 
in the range Eecl < 0.2 GeV is about 32% according to 
the MC simulation. 

The simulation of the Eecl and M^ iss distributions 
in MC is validated using various control samples. The 
Eecl distribution for the B~ — > t~v t signal compo- 
nent is related to beam background and split-off showers 
originating from either E tag or B s - lg decay products. The 
contributions of these components in the signal MC sam- 
ple are 0.04 GeV (16%), 0.12 GcV (50%), and 0.08 GeV 
(34%), respectively. The simulation of the Eecl distri- 



bution is checked with the B~ — > D*°£~D( sample, which 
has a final state similar to the B~ — > t~v t signal if the 
D*° decay products are removed. We also check the dif- 
ference between the detector resolution in data and MC 
for M^ iss with the B~ — > D*°£~I>i sample. We confirm 
that the -Eecl distributions and M^ iss resolutions of data 
and MC are consistent for the B~ — > D*°£~vg sample. 
The background Eecl and M^ iss descriptions by MC are 
checked using sidebands in Mb c and Eecl, events with 
the -Btag reconstructed in a B° mode and events with the 
same B s [ g charge as the B tag . The K\ detection efficiency 
that affects the estimation of the remaining peaking back- 
ground component after the ET£ veto is calibrated using 
the D° —¥ 4>Kg decay sample in data by comparing the 
yields of cj> — > K^Kg and <f> — > K + K~ decays. We con- 
firm the MC expectations for the Eecl and M^ iss shapes 
and verify that the normalization agrees with data after 
the calibrations of the E tag and K\ reconstruction effi- 
ciencies. 

In the final fit, five parameters are allowed to vary: 
the total signal yield and the sum of BB and continuum 
backgrounds for each r~ decay mode. The ratio of the 
BB to the continuum background is fixed to the value 
obtained from MC after the B tag efficiency correction has 
been applied. Other background contributions are fixed 
to the MC expectation. We combine r~ decay modes 
by constraining the ratios of the signal yields to the ra- 
tio of the reconstruction efficiencies obtained from MC 
including the branching fractions of r~ decays (l7j . 

Figure [2] shows the result of the fit to the Eecl and 
M^ iss distributions for all the r~ decay modes combined. 
The signal yield is 62^22 (stat) ± 6(syst), where the first 
and second errors correspond to statistical and system- 
atic uncertainties, respectively. The significance of the 
signal is estimated by y— 2 ln(£o/£max), where £ max 
and Cq are the maximum likelihood and the likelihood 
obtained assuming zero signal yield, respectively. The 
likelihoods are obtained after convolving with a Gaus- 
sian distribution that corresponds to the systematic er- 
ror. We obtain a significance of 3. Oct including system- 
atic uncertainties. The branching fraction is calculated 
by B = N sig /(2eN B+B -), where N sig is the signal yield, 
e is the efficiency, and N B + B - is the number of B + B~ 
events. Equal production of neutral and charged B me- 
son pairs in T(4S) decay is assumed. We obtain 

B(B~ -> t-v t ) = [0.72j£!£(stat) ± O.ll(syst)] x 10" 4 . 

(2) 

The result is summarized in Table HI 

As a check, we fit the Eecl and M^ iss distributions 
floating the yield for each of the four r~ decay modes. 
The resulting yields, as well as the efficiencies and the 
branching fractions are shown in Table HI The branch- 
ing fractions are in good agreement between different t~ 
decays. We also check the result after removing the ET" 
veto, and obtain iV s ig = 65+25 and B(B~ — > t~D t ) = 
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FIG. 2: Distributions of -Eecl (top) and M^ i88 (bottom) 
combined for all the r" decays. The M,^ iss distribution is 
shown for a signal region of -Eecl < 0.2 GeV. The solid 
circles with error bars are data. The solid histograms show 
the projections of the fits. The dashed and dotted histograms 
show the signal and background components, respectively. 



is evaluated by varying the content of each bin by its 
statistical uncertainty. To estimate the systematic error 
due to the possible signal -Eecl shape difference between 
MC and data, the ratio of data to MC for the Eecl his- 
tograms of the B~ — > D*°£~De sample is fitted with a 
first-order polynomial and the signal Eecl PDF is mod- 
ified within the fitted errors. The uncertainties for the 
branching fractions of B decays that peak near zero Eecl 
are estimated by changing the branching fractions in MC 
by their experimental errors [13] if available, or by ±50% 
otherwise. To estimate the uncertainty associated with 
the i?tag efficiency for the signal, B(B~ D*°£~De) ob- 
tained from the B~ — > D*°£~i?( sample is compared to 
the world average value [13] • The results are consistent 
and the uncertainty of the measurement is assigned as 
the systematic error. The uncertainty for the fraction of 
the correctly reconstructed -B tag in the background is ob- 
tained by changing the fractions by errors obtained from 
the Eecl sideband sample. The systematic errors in the 
signal-side efficiencies arise from the uncertainty in track- 
ing efficiency, particle identification efficiency, branch- 
ing fractions of r decays, the reconstruction efficiency 
of 7T°, and MC statistics. The systematic uncertainty 
related to the K £ veto efficiency is estimated from the 
statistical uncertainties of the D° — > 4>Kg control sam- 
ple and the fraction of events with if£ candidates in the 
B~ — > D*°£~vg sample. The estimated systematic errors 
are summarized in Table [TT1 



(0.65t^)xI(T 4 , where errors are statistical only. These 
checks are consistent with the nominal result. In addi- 
tion, we perform one-dimensional fits to Eecl and M^ iss 
and divide the data sample into several subsets. All re- 
sults are in good agreement with the nominal one within 
the statistical errors. 



TABLE I: Results of the fit for signal yields (iV s ig), detection 
efficiencies (e), and branching fractions (B). The efficiencies 
include the branching fractions of the r~ decay modes. The 
errors for iVsig and B are statistical only. 



Decay mode N sig e (1CT 4 ) B (I(T 4 ) 



t — > e u e u T 


16" 


m 

-9 


3.0 


0.68^1? 


T~ — > pTV^Vr 


26' 


i-16 
-14 


3.1 






8j 


10 

8 


1.8 


U - ' -0.59 


T~ — > 1V~1Y V T 


14 


-19 
-16 


3.4 


n 52+ - 72 


Combined 


62 


t-2H 
-22 


11.2 


72+ 27 

u - '^-0.25 



Systematic errors for the measured branching fraction 
are associated with the uncertainties in the signal yield, 
the efficiencies, and the number of B + B~ pairs. The sys- 
tematic error from MC statistics of the PDF histograms 



TABLE 11: Summary of the systematic errors for the branch- 
ing fraction measurement. 



Source 


B syst. error (%) 


Signal PDF 


4.2 


Background PDF 


8.8 


Peaking background 


3.8 


B tag efficiency 


7.1 


Particle identification 


1.0 


7T° efficiency 


0.5 


Tracking efficiency 


0.3 


r branching fraction 


0.6 


MC efficiency statistics 


0.4 


K% efficiency 


7.3 


N B+B - 


1.3 


Total 


14.7 



The branching fraction measured here is lower than the 
previous Belle result with a hadronic tagging method 0] . 
Using the first sample of 449 x 10 6 EE pairs, which corre- 
sponds to the data set used in Ref. @ after reprocessing, 



we obtain B(B 



= [1.08 



+ a35( s t a t)] 



x 10" 



Note that the overlap of events between the two anal- 
yses is small because the reconstruction efficiency has 
increased by more than a factor of three. Assuming 
that all the events used in the previous analysis over- 



7 



lap with present analysis, the remaining events provide 
a result statistically consistent within 1.9er. Using the 
last 323 x 10 6 BB pairs, we obtain B(B~ -> t~D t ) = 
[0.24^0 ^(stat)] x 10~ 4 , which is statistically consistent 
with the result for the first 449 x 10 6 BB data set within 
1.6cr. Our results are also consistent with other publica- 
tions within the errors 

In summary, we measure the branching fraction of the 
decay B~ — > t~v t with hadronic tagging using the fi- 
nal data sample containing 772 x 10 6 BB pairs collected 
by the Belle experiment. We find evidence of B~ — > 
t~v t with a branching fraction of [0.72to'2s( s t a t) ± 
O.ll(syst)] x 1(T 4 , and a signal significance of 3. Oct in- 
cluding systematic uncertainties. By employing a neural 
network-based method for hadronic tagging and a two- 
dimensional fit for signal extraction, along with a larger 
data sample, both statistical and systematic precisions 
are significantly improved compared to the previous anal- 
ysis [6j. The result presented in this paper supersedes 
the previous result reported in Ref. @. Combined with 
the Belle measurement based on a semileptonic B tag- 
ging method Q, the branching fraction is found to be 
B(B~ -> T~v r ) = [0.96±0.22(stat)±0.13(syst)] x 10" 4 , 
with a 4. Oct signal significance including systematic un- 
certainties. The measured branching fraction is now con- 
sistent with the SM expectation obtained from other ex- 
perimental constraints. Using this result and parameters 
found in Ref. [13], we obtain f B \V ub \ = [7.4 ± 0.8(stat) ± 
0.5(syst)] x 10~ 4 GcV. Our result provides stringent 
constraints on various models including charged Higgs 
bosons. 
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